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We present globally regular vortex-type solutions for a pure SU(2) Yang-Mills 
field coupled to gravity in 3+1 dimensions. These gravitating vortices are static, 
cylindrically symmetric and purely magnetic, and they support a non-zero chromo- 
magnetic flux through their cross section. In addition, they carry a constant non- 



Abelian current, and so in some sense they are analogs of the superconducting cosmic 
strings. They have a compact central core dominated by a longitudinal magnetic 
field and endowed with an approximately Melvin geometry. This magnetic field 

o 

component gets color screened in the exterior part of the core, outside of which the 



fields approach exponentially fast those of the electrovacuum Bonnor solutions with 
a circular magnetic field. In the far field zone the solutions are not asymptotically 
flat but tend to vacuum Kasner metrics. 

PACS numbers: 04.20.-q, 04.25.Dm, ll.15.-q, 11.27.+d 

Classical solutions in non-Abelian gauge models coupled to gravity in four space-time 
dimensions attracted a lot of attention after the discovery by Bartnik and McKinnon 
of particle-like solutions in the Einstein- Yang-Mills (EYM) theory and after the construc- 
tion of the corresponding black holes (see [2| for a review). Since then a variety of static 
and stationary axisymmetric EYM solutions have been considered, but static, cylindrically 
symmetric EYM configurations have remained unexplored for some reason. To fill this gap, 
we consider in this letter static, cylindrically symmetric, purely magnetic, globally regular 
four- dimensional EYM solutions of cosmic string type. 

Surprisingly, we find strings endowed with a constant current along their symmetry axis, 
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and so in some sense they can be considered as analogs of Witten's superconducting cosmic 
strings jjj. In the latter case the superconductivity arises due to the presence of two complex 
scalars in the theory In our case their role is played by two color components of the self- 
interacting Yang-Mills field, which behave exactly in the same way as the scalars in Witten's 
model. Specifically, one of them develops a non-zero condensate value in the vortex core 
and vanishes at infinity, while the other one behaves the other way round. Since the gauge 
symmetry of our theory is not broken, the gauge field coupled to the current is long-range. 
The slow fall-off of the energy density in the direction orthogonal to the string then does 
not let the solutions to be asymptotically flat, although they are asymptotically Ricci flat. 
This is the principal difference of the EYM strings as compared to the previously studied 
currentless Abelian self-gravitating cosmic strings of the Nielsen-Olesen type 

In what follows we present a numerical evidence for the existence of a family of globally 
regular EYM string solutions labeled by an integer winding number v and by a real parameter 
p determining the value of the magnetic field at the symmetry axis. These EYM strings 
are filled with a self-interacting chromo-magnetic field giving rise to a constant current 
along the string symmetry axis and also to a chromo-magnetic flux through the string 
cross section. The flux and current can be defined in a gauge-invariant way within the 
approach of asymptotic symmetries jfj . The geometry of the solutions interpolates between 
flat Minkowskian geometry at the symmetry axis and vacuum Kasner geometry in the far 
field zone. Interestingly, these solutions survive even if the coupling to gravity is off and 
the metric is flat everywhere. They describe then 'pure Yang-Mills superconducting strings' 
with a regular central core outside of which there remains only a U(l) component of the 
gauge field, which, being coupled to the current, diverges logarithmically for large r. As 
a result, the total energy per unit length for the flat space solutions diverges at large r. 
Getting back to curved space, this divergence is cured by gravity. 

We consider the SU(2) Einstein- Yang-Mills (EYM) theory 

where F" = d^A^ — d v A a + e a b c A h A c u and e is the gauge coupling constant. With A = 
^TaA^dx^ where r a are the Pauli matrices, the local SU(2) gauge symmetry of the action 
is expressed by A — > U(A + id)U~ l where U(x) G SU(2). We shall be interested in static, 
cylindrically symmetric systems invariant under the action of the three commuting Killing 
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vectors d t , d^, d z . The spacetime metric then can be parameterized as 

ds 2 = l 2 {N(r) 2 dt 2 - H(r) 2 dr 2 - L(r) 2 d<p 2 - K(r) 2 dz 2 }, (2) 

where I is the length scale and all the other quantities are dimensionless. For the gauge field 
we make the following purely magnetic ansatz, 

A = r 2 R(r) dz + r 3 P(r) dip. (3) 

Using Eqs.Q,© we obtain S = J d 4 xC where k = 8nG/(el) 2 and 
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We notice that the reduced Lagrangian C is invariant under P — > aP, L — > aL, R — > PR, 
K — > (3K, N — > N ', H — > H , with constant a, [3. It also admits a discrete symmetry P «-> R, 
L <r-> K, which is equivalent to z <-> (p. Varying £ gives the field equations, which can be 
represented in the form 

fN' 2 L'S' L' 2 P 2 R 2 \ ^ 
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plus a first order constraint equation 



N'S' L'S' L' 2 _ k fP' 2 L 2 R' 2 N 2 P 2 R 2 * 
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Here S = LK and, after varying, we have imposed the coordinate condition H = N. The 
energy density is 

2N 2 L 2 2N 2 K 2 2S 2 ' 1 ' 

while the other non-trivial components of T£ are obtained by choosing the signs in front of 
the three terms in this formula as T r r ~ (— , — , +), T 2 Z ~ (+, — , — ), ~ (— , +, — ), such 
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that Tjf = 0. Associated with the two continuous global symmetries of £ there are two 
Noether charges 

(N> V PP'\ 
„(N' S' V L 2 RR'\ 

whose conservation can be checked straightforwardly. 

We want the fields to be regular at the symmetry axis r = 0. This implies that one 
should have P(0) = v E Z, since in this case, with U = e - ^ 1 " 3 ^, one can pass to the gauge 
where the gauge field will be regular at the axis: 2A = Ut 2 \J~ x Rdz + t% (P — v) dip. The 
metric should become Minkowskian for r — ■> 0. The most general local power series solution 
of Eqs.(jnj) with such boundary conditions reads 

P = v-pr 2 + 0(r 4 ), R = qr v + 0{r u+2 ), (9) 
N = l + 0(r 2 ), S = r + 0(r 3 ), L = r + 0(r 3 ), 

where p, q are two integration constants. This solution fulfills also the constraint ( J^j l. In- 
serting (jUJ) to (JBJ) fixes the values of the Noether charges: Qi = 2vnp — 1 and Q 2 = 0. We 
can now integrate the equations starting from the origin on, with the boundary condition 
©, and this reveals the following picture. Fixing a value of p and varying q, the numerical 
solution generically terminates at a finite point r = r* (q) where P' becomes large and either 
positive or negative, depending on q. Adjusting properly the value of q, we can postpone 
further and further the moment r*(g) where the numerical procedure crashes. We thus ex- 
tend the solution from the axis farther and farther to the asymptotic region, and we observe 
then that the amplitude P approaches zero very quickly with growing r. As a result, we 
conclude that at large r the globally regular solutions approach configurations with P = 0. 

Setting P = in Eqs.(jHJ), the SU(2) gauge field configuration becomes U(l), in which 
case the general solution of Eqs.(jHJ, let us call it N , L , S , R , reads (f| 

L = ai x 1+m + a 2 x x ~ m , N = a 3 x m2 -%, (10) 
S = a 4 x, R = a 5 + a 6 — — , x = r — r , 

where m, a 1; . . . a 6 , r are 8 integration constants. Imposing the constraint © would imply 
also that a 6v /Ka7 = ±a 4V /2o2, but we do not demand this at the moment. We conclude that 
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FIG. 1: The globally regular solution of Eqs.(JSJ) with p = 0.8, v = 1. The central core extends 
up to £ rj 1.2, after which the configuration practically coincides with the Abelian solution (jlUj) of 
Bonnor with m = —1.94. For £ > 3 the gauge field dies out and the solution becomes Kasner. 

at large r our solution is given by 

L = L + 8L, N = N + 8N, S = S + 5S, 

R = R + 5R, P = 5P 1 (11) 

where the deviations 5L, . . . , 5P vanish as r — >• oo. To determine them we insert (jllj) to 
Eqs.flSJ) and linearize with respect to the deviations. This gives 

xr> ( L N R \ . 

7 7^ exp v — — 7 

where aj is a new integration constant and the dots stand for the subleading terms. It follows 
also that SN, 5L, 5S, and 5R are all of the order {SP) 2 and no new integration constants 
appear in these variations, since the background solution N , L , S , Ro already contains the 
maximal number of integration constants. 

Summarizing, the local solution at large r contains 9 free parameters and is given by 
Eqs. (llO|) - (fT2|) . while that at small r contains 2 integration constants in Eq.©. We then 
integrate numerically Eqs.(jHJ) to extend these local solutions to the intermediate region where 
we match them. This gives us global solutions in the interval r e [0, oo). To fulfill the 10 
matching conditions for the 5 functions in Eqs.(j5J) and for their first derivatives, we have 
in our disposal 2+9=11 free parameters, and so there is one free parameter left after the 
matching. This parameter, we choose it to be p in Eq.(|9~|). labels the global solutions. 

No special care is to be taken about the constraint (JBJ), since being imposed at r = 
it 'propagates' to all values of r. As a result, the values of the asymptotic parameters of 
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FIG. 2: The parameter m{v,p), the Tolman mass M(v,p) and the current I(is,p) = I^/ftir) for 
the globally regular solutions with v = 1,2. 



the global solutions found by the matching will automatically satisfy a^^nai = ±a4v2«2- 
Similarly, the axis values of the Noether charges Qi = 2unp — 1, Q 2 = propagate to the 
large r region and impose the asymptotic conditions 



a^[m 2 — 1) = Qi, (2 + m)a 4 + 2^2na 1 a 2 a 5 = 0. 



(13) 



We thus obtain a family of globally regular solutions of Eqs.flni) labeled by a real parameter 
p > and by an integer z/, and we use the length scale for which k = 2. The typical solution 
is shown in Fig.l. These EYM strings have a compact central core filled with a chromo- 
magnetic field, outside of which P pa and the remaining gauge field becomes effectively 
Abelian. The interior par of the core, where R ~ 0, is also effectively Abelian and can be 
well described by the Melvin solution. The latter can be obtained from the general solution 
(jl(J|) by applying the z (p symmetry and adjusting the constants such that the boundary 
conditions at r = are fulfilled: 

k o o T pr 2 , , 

N = K = l + -p 2 r 2 , L = -, P = u- P —. (14) 

Since this inner solution has a magnetic field along the z axis, we call it z-string. The metric 
in this limit acquires an additional symmetry with respect to Lorentz boosts in the (0, z) 
plane, and the energy- momentum tensor in this limit is such that Tq = T 2 2 = — T r r = — Tg. 

The outer part of the central core is the transition region where the non-Abelian effects 
become essential, since the R field starts growing there, playing a role of the effective mass 
term for the P field, and so the latter starts falling down to zero exponentially fast. As a 
result, outside of the central core one has P — » and the whole configuration approaches 
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FIG. 3: The Euclidean 2-geometry of the (r, <p) string orthogonal section can be realized on a 2- 
surface of revolution in three dimensional Euclidean space obtained by rotating the shown planar 
curve around the Z-axis. The curve shown corresponds to the string with v = l,p = 0.75. 

exponentially fast one of the Abelian Bonnor solutions ([10)1 . These, in their turn, contain 
an intermediate region filled with a circular magnetic field, which we call 0-string and which 
has Tq = T£ = — T r r = — T 2 Z . This magnetic field falls off polynomially for large r, and in the 
far field zone the solutions approach asymptotically vacuum Kasner metrics. Summarizing, 
the whole solution can be viewed as a non-linear superposition of a z-string placed at small 
r and a 0-string located at larger r, which is demonstrated by the profile of the radial energy 
density, SN 2 Tq, clearly showing two peaks (see Fig.l). 

It is natural to wonder if the solutions can be asymptotically flat, as in the case of 
gravitating Nielsen-Olesen strings 4]. However, taking a linear combination of the first 
equation in (j5J) with the constraint (jUJ) and integrating gives 



The integral on the right here is manifestly positive and equals the Tolman mass, M, while 
the left hand side is equal to (m 2 + |m|)a 4 . As a result, one always has m ^ and so 
solutions with Tq 7^ cannot approach flat metric for large r. 

Apart from the mass M, another important parameter of the EYM strings is the current 
flowing along them. Although there is no Higgs field in the theory, it is still possible to in- 
troduce conserved and gauge invariant currents within the asymptotic symmetries approach 
of Abbott and Deser [5j. Specifically, since our solutions are asymptotically vacuum, one 
has A = A°° + a, where A°° = iUoodll^ 1 and a — >• as r — > 00, but otherwise a need not to 




(15) 
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be small. Defining T>^a v = V fl a i 



i[A°f,a v ] one can introduce 



(16) 



and then the conserved and gauge invariant currents are J£ = V a f^ ,a . The total current 
through the string cross section is then given by 



Our strings thus carry a persistent current, and for this reason they can be called 'super- 
conducting'. 

At the axis we have similarly A = A + a, where A = iUqcLUq 1 and a — > for r — > 0. 
Defining then //£ as in (jlfij) but with 'oo' replaced by '0', we can calculate its flux through 
the sting cross section, S, 



and this is 'quantized' as for the Nielsen- Olesen vortex. 

The mass, current and the parameter m as functions of p are shown in Fig. 2 for solutions 
with v = 1,2. We see that m is always negative and tends to zero for p — * 0, in which limit 
the gauge field becomes pure gauge with R(r) = and P(r) = u, so its energy vanishes 
and the metric becomes flat. For p = 1/(2kv), when Qi = 0, Eq. (jl3|) implies that m = — 1, 
in which case #oo becomes proportional to g w for large r, and so the metric acquires the 
asymptotical boost symmetry in the (t, (p) plane. 

Another interesting limit is k — > with fixed p, in which case the gravity switches off and 
the geometry becomes flat. The gauge field, however, remains non-trivial. Such Minkowski 
space solutions also show a regular central core where R ~ and P^O, outside of which P 
is exponentially small, but the large r behavior of R is now different. Specifically, for small 
k one has m = -^/k/2I, ai — 1/2 = 0(y/n), <22 — >• 1 — cti, using which in Eq. (|10j) and 
taking the limit gives R = Cq +Xln(r) for large r. Of course, the same results are obtained 
by solving the two Yang-Mills equations in (JHJ) with N = 1, L = S = r. The amplitude R 
is thus now logarithmically divergent at infinity, as it should be, since we effectively have 
here a U(l) gauge field coupled to the conserved current. We thus obtain 'superconducting 
strings' made of pure Yang-Mills field in flat space, and these, to the best of our knowledge, 
have also never been described in the literature. Getting back to the self-gravitating case, 




(17) 




(18) 
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one can say that the gravity provides a 'dynamical cut-off' by rendering R finite at large r, 
which makes the total energy finite. 

It is interesting to visualize the two-dimensional geometry dl 2 = N 2 dr 2 + L 2 d(p 2 of the 
(r, tp) plane orthogonal to the string in terms of embeddings. The same geometry can be 
realized on a 2-surface of revolution in three dimensional Euclidean space with coordinates 
X, Y, Z obtained by rotating around the Z-axis the planar curve defined parametrically by 
relations X = ±L(r), Z = Z(r) with Z' 2 = N 2 — L' 2 . This defines a 'flowerpot' surface 
shown in Fig. 3. We observe that the circumference of a circle centered at the string axis in 
the plane orthogonal to the string grows up slower than the radius of the circle. 

Summarizing, we have presented superconducting cosmic string solutions for a pure self- 
gravitating Yang-Mills field. We have not studied their stability yet. As a matter of fact, 
many solutions in the EYM theory are unstable. However, since we now have the conserved 
current, there is a chance that it may stabilize the solutions. Other interesting problems 
to study would be to consider the black hole generalizations for our solutions or to include 
other fields. We also notice that replacing in all above formulas t <-> z gives solutions with 
an electric field. They do not have a current, but possess instead a chromo-electric charge. 
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